SINGULAR PERTURBATIONS OF MEAN CURVATURE FLOW 

GIOVANNI BELLETTINI, CARLO MANTEGAZZA, AND MATTEO NOVAGA 

Abstract. We introduce a regularization method for mean curvature flow of a subman- 
ifold of arbitrary codimension in the Euclidean space, through higher order equations. 
We prove that the regularized problems converge to the mean curvature flow for all times 
before the first singularity. 
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1. Introduction 

It is well known that a smooth compact submanifold of the Euclidean space, flowing 
by mean curvature, develops singularities in finite time. This is a common aspect of 
geometric evolutions, and motivates the study of the flow past singularities. Concerning 
the mean curvature motion, several notions of weak solutions have been proposed, after the 
pioneering work of Brakke , see for instance [H El El El El UHl HU [T31 1201 12U 1221 12^11221 
We recall that some of these solutions may differ, in particular in presence of the so-called 
fattening phenomenon (see for instance 

Following a suggestion of De Giorgi in [T2j, we introduce and study a regularization of 
mean curvature flow with a singular perturbation of higher order, which could lead to a 
new definition of generalized solution in any dimension and codimension. 

Let us state our main result. 
Let if : M —>■ M n+m be a smooth compact n-dimensional immersion in IR n+m . For k > 
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[n/2] +2 (where [n/2] denotes the integer part of n/2) and e > we consider the functional 

/ {l+e\A k \ 2 ) dfi, 

J M 

where \x is the canonical volume measure associate with the metric g induced on M via the 
immersion (p. With A k we denote the ^-differential in M. n+m of the function A M given by 

A M (xj = l x l 2 ~ [d M (x)] 2 

where [d M ] 2 is the square of the distance function from ip(M), which is smooth in a neigh- 
borhood of a point of the submanifold without self-intersections. Since locally on M every 
immersion is an embedding, we can define A k also at such points. More precisely, the 
tensor A k is defined as 

A u d k A M 



dx n ...dx lk 



for every fc-uple of indexes i±, . . . , z& G {1, . . . , n + m}. 

We remark (see |T4^ Prop. 2.2 and Cor. 2.4]) that the tensors A k and V fc ~ 3 B, where B is 
the second fundamental form of M and V is the covariant derivative associated with the 
induced metric g, are strictly related, hence, in a way the functional Q\ is a perturbation 
of the area functional by a term containing the squares of the high order derivatives of the 
curvatures of M. 

c By means of Theorem 4.5 and Theorem 5.9 in j2] and the results of [Tl], the gradient 
flow associated with the functional Q e k is given by the PDE system 

(k — 2)-times 

(1.1) ^- =H + 2 £ fc(-lf(A M oA M o ... o A M Hj +eLOT 

where H is the mean curvature vector and LOT denotes terms of lower order in the curva- 
ture and its derivatives. 

We can see then that (|l.lj) is a singular perturbation of the mean curvature flow, and co- 
incides with it when e = 0. In (see also jTS] and j2Hl) it is proved that for every e > 
the system in (|1.1|) admits a unique smooth solution defined for all times; we are then 
interested in the convergence to the mean curvature flow when e — > 0. 
Our main result is the following. 

Theorem 1.1. Let (p : M — > ]R n+m be a smooth immersion of a compact n-dimensional 
manifold without boundary. Let T sing > be the first singularity time of the mean curvature 
flow ip : M x [0,T sing ) -> W n+m of M. For any e > let <p £ : M x [0, +oo) -> R n+m be 
the flows associated with the functionals Q e k , with k > [n/2] + 2, all starting from the same 
initial immersion ipo- Then the maps (f £ converge locally in C°°(M x [0,T S i ng )) to the map 
ip, as e — > 0. 

Example 1.2. In case of immersed plane curves 7 : S 1 — > IR 2 (n = m = 1) the simplest 
choice is k = 3. Since it turns out that |^4 3 | 2 = 3k 2 , where k is the curvature of 7, in this 
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simple case the approximating functionals read as 



(l + en 2 ) ds 



7 



where s is the arclength parameter, and we have replaced 3e with e. The regularized 
system which approximate the curve shortening flow is then 



where v is a suitable choice of the normal unit vector to the curve. 

The crucial point in order to prove Theorem 1 1.1 1 is to obtain e-independent estimates of 
the curvature and its derivatives in order to gain sufficient compactness properties. We get 
these by computing the evolution equations satisfied by the L 2 norms of the derivatives of 
the second fundamental form of the flowing manifolds, and by estimating via Gagliardo- 
Nirenberg interpolation inequalities. 

At the present moment we are not able to characterize the limit of the approximating flows 
after the first singularity, as the proof of Theorem 11.11 relies heavily on the smoothness of 
the mean curvature flow in the time interval of existence. Our goal would be to provide 
some limit flow defined for all times, thus providing a new weak definition of solution in 
any dimension and codimension. 

We mention the simplest open problem in defining a limit flow after the first singularity. 
It is well known (Gage-Hamilton ^HJ EI and Huisken JH]) that a convex curve in the 
plane (or hypersurface in M ra+1 ) moving by mean curvature shrinks to a point in finite 
time, becoming exponentially round. In this case we expect that the approximating flows 
converge (in a way to be made precise) to such point for every time after the extinction 
one. 

The plan of the paper is the following. In Section |2] we give some notation and we recall 
the relations between the squared distance function and the second fundamental form and 
its covariant derivatives. In Section El in order to make the line of proof clearer, we work 
out in detail the e-independent estimates in the simplest case of plane immersed curves; 
also in this special case, the result appears to be nontrivial. In Sections |H El and El we 
consider the general case of a n-dimensional submanifold of M ra+m . Section [7| is devoted to 
show Theorem ll.il 

We remark here (but we will not discuss such an extension in this paper) that our method 
works in general for any geometric evolution of submanifolds in a Riemannian manifold 
till the first singularity time, even when the equations are of high order (like, for instance, 
in the Willmore flow, see [23 EH EH] ) , choosing a regularizing term of appropriately high 
order. 

Finally, it should be noted that, looking at the evolution equation (jl.lj) . these pertur- 
bations of the mean curvature flow could be considered, in the framework of geometric 
evolution problems, as an analogue of the so-called vanishing viscosity method. Indeed, 
we perturb the mean curvature flow equation with a regularizing higher order term multi- 
plied by a small parameter e > 0. The lower order terms, denoted by LOT, which appear 




$7 



(k — 2ed 2 n — en 3 ) v. 
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in (jl.lj) are due to the fact that we actually perturb the area functional and not directly 
the evolution equation. However, the analogy with the classical viscosity method cannot 
be pushed too far. For instance, because of the condition k > [n/2] + 2, our regularized 
equations are of order not less than four (precisely at least four for evolving curves, at 
least six for evolving surfaces). Moreover, as the Laplacians appearing in equation 
are relative to the induced metric, the system is quasilinear and the lower order terms are 
non linear (polynomial). 



We denote with e\, . . . ,e n+m the canonical basis of W a+m and with ( , ) its standard 
scalar product. 

We let M C M n+m be a smooth, compact, n-dimensional, regular submanifold without 
boundary, then T X M, N X M C M n+m are, respectively, the tangent space and the normal 
space to M at x G M C R n+m . 

The distance function d M and the squared distance function n M from M are given by 



for any x G R n+m (we will drop the superscript M when no ambiguity is possible). In this 
section we recall some facts from [2. and about the distance function and the relations 
between the high derivatives of r] M and the second fundamental form of M. 

When M is embedded, there exists an open neighborhood Q C M. n+m of M such that 
d is smooth in Q \ M and rj M is smooth in the whole of Q. If M is only immersed, at any 
point of M we consider the distance (we still use the symbols d M and rj M for simplicity) 
from an embedded image of a suitable neighborhood of the point; in this case the regularity 
properties of d M and r\ M hold in a neighborhood (still denoted by Q) of such an embedded 
image. 

Clearly, i] M and V?7 M (x) = at every x G M; moreover, for every x G f2 we have that 
x — Vr) M (x)/2 is the unique point in M of minimum distance from x (the projection of x 
on M), that we denote with tt m (x). 

Another nice property of the squared distance is that, for every x G M the Hessian matrix 
V 2 rj M (x) is twice the matrix of orthogonal projection onto the normal space N X M. We 
will denote respectively with X M and X 1 - the orthogonal projections of a vector X on the 
tangent and normal space of M. 

Let x G M and 1,7 6 T X M, the vector valued second fundamental form of M at the 
point x is given by 



2. Notation and Preliminaries 



x) = inf \x — 



and 



V M ( X ) = [d M {x)f 




where we extended locally the two vectors X, Y to tangent vector fields on M (the deriv- 
ative is well defined since X is a tangent vector at x). 
If { z/ a}o=i,...,m is a local basis of the normal bundle we have 
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where here and throughout all the paper we use the convention of summing on repeated 
indices. 

We will see B as a bilinear map from T X M x T X M to M n+,n , hence, as a family of n + m 
bilinear forms B fc = (B, e^) : T X M x T X M — > R. Moreover, we consider B acting also on 
vectors of R n+m , not necessarily tangent, by setting B(V,W) = B(V M ,W M ) G N X M C 
R n+ m for every pair v,We R n+m . With such definition, B£- = (B(e i; e 3 -), e k ). 
It is well known that B is a symmetric bilinear form and its trace is the mean curvature of 
components H fe = B k 



3d 

We introduce now the function 



A M (x) 



\x\ 2 - [d M {x)} 2 



2 

smooth as i] M in the neighborhood Q of M, and we set 

a m (x) 

for the derivatives of A M at every point x G Q. 

The following Proposition (see [2] for the proof) shows the first connection between the 
second fundamental form and the function A M (or equivalently, the squared distance func- 
tion). 

Proposition 2.1. The following relations hold. 

• For any x G Q, the point VA M (x) is the projection point ir M (x). 

• If x G M, then V 2 A M (x) is the matrix of orthogonal projection on T X M. 

• For every x G M, 

r>k _ aM (r A M \ 

Afjk — + B* fc + B J ki , 



We define now the fc-derivative tensor A k (x) acting on fc-uples of vectors Xi G TOn+m 



■ 



where Xj = X\e^ as follows 

A h (x)(X 1 ,...,X h ) = AZ_ jik (x)X?...Xj- i 
notice that the tensors A k are symmetric. 

For notational simplicity, we drop the superscript M on A k ; for the same reason, we also 
avoid to indicate the point x G M in the sequel. 

The tensors A k and V fe ~ 3 B are strictly related by a recurrence formula proved in [TH 
Prop. 2.2 and Cor. 2.4]. 

Remark 2.2. We underline here an important convention used in the paper. Due to the 
high codimension, we will work with several tensors (like normal vector fields or the second 
fundamental form B) taking values in IR n+m ; these tensors will be considered as families of 
n + m tensors with values in M. With this convention, for instance, VB means that we are 
considering the family of covariant derivatives of the tensors B 1 , . . . , B n+m , one component 
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at time. Unless otherwise specified, this convention will be used even also for tangent 
vector fields, that is, when X is a tangent vector field, VX is not the covariant derivative 
of X but the derivative of its components in the basis of M n+m . 

In all the paper we write T * S, following Hamilton [1$; , to denote a tensor formed by 
contraction on some indices of the tensors T and S using the coefficients g % K 
If T\, ... ,Ti are tensors (here I is not an index of the tensor T), with the symbol 

i 

® Ti 
i=i 

we mean T\ * T% * • • • * 7} . 

Definition 2.3. We use the symbol p s (V'B) for a polynomial (with the * product) tensor 
with constant coefficients in the coordinate basis d<p/dxi, the second fundamental form B 
and its derivatives up to the order I at most, such that every of its monomial is of the form 

N 

® V jk B with < j k < I and iV > 1 

k=l 

or 

N 

® V jk E . 

k=l ax 

where, in both cases, the rescaling order s equals 

N 
k=l 

We use instead the symbol q s (V'B) for a polynomial of the kind p s (V'B) such that the 
contraction with the metric is total, both in the covariant and in the M n+m -indices. 
As the contraction in the ambient space R n+m "cancels" all the basis elements dip/dxi 
appearing in the formulae, it follows that every monomial of q s (V'B) has the form 



V ifc B-^ with < j k < I and iV > 1 



N N 

© 
k=l 



V ifc B with < j k < I and ^(jk + 1) = s, 



k=i 

where the covariant indexes are all completely contracted with g % K 

Remark 2.4. See the paper [2E1 Sect. 2] for more details on these polynomials and the 
geometric interpretation of the rescaling order. Notice that, differently from J2H|, here we 
need to consider in p s (V'B) monomials of two types, because of the codimension higher 
than one. 



We advise the reader that the polynomials p s and q s may vary from line to line, and 
similarly the constants (usually indicated by C). 
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3. Evolving Plane Curves 

Let 7 G C°°(S l ] M 2 ) be a regular immersed closed curve in the plane M. 2 . Let r = 
lx/\lx\ = 7s and v = Rr be respectively the tangent and the normal to the curve 7, where 
R is the counterclockwise rotation of tt/2 in the plane, and 7^ = d x ^y. 
We recall that d s = d x /\j x \ and 

(3.1) d s r = ku, d s v = —kt 

where k is the curvature of 7. In the sequel we let L = L(7) = J Ids be the length of the 
curve 7. 

Let us consider the functional 



Q £ {l) = I (I+ek 2 ) ds, 

J ■y 



which is obtained from Q\ (with n = m = 1) by replacing 3e with e. Set 

W = -k + 2ed 2 s n + en 3 . 

Then the gradient flow by Q £ is given by a smooth map 7 : S 1 x [0, +00) — > R 2 which is an 
immersion for any t G [0, +00), equals a given immersion 70 at time t — 0, and satisfies 

(3.2) d a = -Wu, 

where d t — J^. For notational simplicity, we omit the dependence of 7 on e. 
Lemma 3.1. We have 

d s d a = -(d s E £ )lS + KE £ T, 

in particular 

(3.3) (d s d t J, r) = kE £ . 

Proof. It follows from equations (|3.1|) and the evolution equation ()3.2|) . □ 
Lemma 3.2. Let ^ be a smooth closed curve, then 

(3.4) l^^J/ ds - 

Proof. By Borsuk and Schwartz-Holder inequalities we have 

1/2 

2vr < / \k\ ds < i I n 2 ds ) L 1/2 . 

□ 



7 



Lemma 3.3. The following commutation rule holds: 
(3.5) dtd a = d s d t - KE £ d s . 
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Proof. Observing that j^j = j^jdt = d s df, we have 

dd —d ( — ^ - dtdx - ® tlxl ® x - —d - I 7x ® tlx ^ dx 



I7x|/ 1 7x| I7x| d I7x| \|7x| 17x1/ |7x| 

= d s d t - (r,d s da)d s . 

Then the commutation rule ()3.5j) follows from equation ()3.3j) . □ 
Lemma 3.4. We have 

(3.6) d t n = -d 2 E £ - k 2 E £ = d 2 s K + k 3 - 2e3 a s k - QsK(d s K) 2 - 5eK 2 d 2 K - en 5 . 

Proof. We have 

d t K = d t (d s r, v) = (d t d s r, v) . 
Therefore, using formula ()3.5j) we have 

d t K = {d s d t d sll U) ~ KE £ (d s T, U) = (d 2 d t l, U) - (d s [KE £ d sl ] } U) ~ K 2 E £ . 

Using the evolution law (|3.2|) we get 

(d 2 d a , v) = -(d 2 {E £ u), v) = -d 2 E £ + E £ (d s {KT), v) = -d 2 E £ + k 2 E £ . 
In addition, 

(d s [KE £ d s -f], v) = KE £ (d s r, v) = k 2 E £ . 
Hence dtK = —d 2 s E e — k 2 E £ and the last equality in ()3.6j) follows by expanding E e . □ 



k\ 



Remark 3.5. For e = 0, formula (|3.6|) gives the well known evolution equation K t = d 2 K+ 
valid for motion by curvature, see fTJ Lemma 3.1.6]. 

We recall now the following interpolation inequalities for closed curves, see pag. 93]. 

Proposition 3.6. Let 7 be a regular closed curve in M 2 with finite length L. Let u be a 
smooth function defined on 7, m > 1 and p G [2, +00]. If n E {0, . . . , m — 1} we have the 
estimates 

(3.7) \\d: u \\ LP < c^jd^uiuniT + ^h\\ L , , 

where 

n + I/2-I/0 . s 

a = '- f -t e 0, 1 

m 

and the constants C n ^ p and B n ^ p are independent of '7. 

Clearly inequalities ()3.7|) hold with uniform constants if applied to a family of curves 
having lengths uniformly bounded below by some positive value. 

Remark 3.7. In the special case p = +00, we have a = n+1 ' 2 ; and 



ion II ^ /~t II am lie II ill— c 1 B nm ., 
\d s u\\ L ™ < C n,mWs U\\ L2 \\U\\ L2 + J-^\\U\\ L2 
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Remark 3.8. In the particular case n = 0, m = 2, p = 6 we get a = 1/6 and 

2 I 5 C 

\\u\\ L6 < C\\d s u\\hW u Wh + ttIM| L 2 > 

L3 

for some C > 0, hence, by means of Young inequality \xy\ < ^\x\ a + l\y\ b , 1/a + 1/6 



choosing a = 6 = 2, x = v^H^-uH^, 2 and y = -^||u||^3, we obtain 



(3.8) J^u 6 ds< J{d 2 s u) 2 ds + C (^J^u 2 ds^j + 9_(^Ju 2 d^j 

In the particular case n = 0, m = 1, p = 4 we get cr = 1/4 and 

I ..a C 
\\u\\ L 4 < C\\d s u\\l 2 \\u\\ 4 L2 + -t||m|| L 2 , 

hi 

hence, reasoning as before, 

(3.9) J^u 4 ds< J(d s u) 2 ds + c(^J^u 2 ds^j +^Q^ 2rfs ) • 

We are now ready for the estimates. We recall that 

(3.10) d t ds = kE £ ds = (-k 2 + 2ek8 2 k + ek 4 ) ds. 
Lemma 3.9. We have 

(3.11) d t f K 2 ds= I ( - 2{8 s k) 2 + k 4 - Ae(d 2 K) - ek 6 - Aek 3 8 2 k) ds . 

J-y J *y 

Proof. From ()3.1()j) and Lemma [3.41 we get 

d t / n 2 ds=2 / Kd t nds+ (-k 4 + 2ek 3 3 2 k + ek 6 ) ds 

J 7 J 7 J 'y 

= 2 k{8 2 s k + k 3 - 2e0 4 k - Qek(3 s k) 2 - hEK 2 d 2 s n - EK b ) ds 

J-y 

+ / (-K 4 + 2EK 3 8 2 S K + EK & ) ds 
J-y 

= / {2Kd 2 s K + K 4 -AEKd 4 s K-l2EK 2 (d s K) 2 -%EK 3 d 2 s K-EK & )ds. 

Therefore, integrating by parts, we obtain 

d t f n 2 ds= I ( - 2(d s n) 2 + k 4 - 4:E(d 2 n) 2 - ek 6 - 12ek 2 (0 s k) 2 - 8EK 3 d 2 K) ds 

= I ( - 2(d s n) 2 + k 4 - 4:E(d 2 n) 2 - ek 6 - 4:EK 3 d 2 n) ds, 

J-y 

where in the last equality we used the fact that —3 J K 2 (d s n) 2 ds = J n 3 d 2 nds. 
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Proposition 3.10. The following estimate holds 

(3.12) 8 t J^ 2 ds<C (^J^ 2 ds^j + c(^J^ 2 ds^j , 

where C is a constant independent of e. 

Proof. Adding to the right hand side of equation (|3.11|) the positive quantity 2e(d 2 n + k 3 ) 2 
we get 

8 t f K 2 ds< j ( - 2{d s K) 2 + k 4 - 2e(d 2 K) 2 + en 6 ) ds . 
Using now inequalities (|3.8|) and (j3.9j) we obtain 

d t J k 2 < J ( - {8 s k) 2 - e[d 2 s K) 2 ) ds + Ce (j k 2 ds^j +jj(J^ 2 ^ 

<C^j^K 2 ds^j +c(^J^ 2 ds^j +9-Q K 2 ds^ +^( ln 2 ds 
<c( [ n 2 ds] + c( [ K 2 ds^ 



where we supposed e < 1 and in the last inequality we used the geometric estimate (|3.4j) . 

□ 

We deal now with the higher derivatives of the curvature. 
Since here we are working in dimension and codimension one, for the rest of this section all 
polynomials in the curvature k and its derivatives are completely contracted, that is they 
belong to the "family" q r (d l s n) (see Definition I2.3j) : moreover, every of their monomials is 
of the form 

N 

Y[ 0*« with < j, < I and N > 1 

8=1 

with 

JV 

r 

i=i 

as the * product in this case is simply the usual product. 



Lemma 3.11. For any j G N the following formula holds: 

(3.13) dtd> a K = di +2 K + q i+3 (d>) - 2ed{ +4 K - 5en 2 di +2 n + ec\ j+ \d j s +1 k) . 
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Proof. We argue by induction on j. 

The case j = in ()3.13|) is equation Ij3.6|) . where q 5 (d s n) = — 6/t(<9 s /t) 2 — k 5 . 
Suppose that (J3.13|) holds for (j — 1); using the commutation rule ()3.5j) we get 

9fc£^'/c = d s d t d J s ~ 1 K + k(k — 2ed 2 s K — £k z )&{k 

= d s [di +1 K + q 3 " +2 (52'- 1 «) - 2 £ df' +3 K - 5e« 2 a£ +1 K + eq 3 ' +4 (52«)] 
+ q^ +3 (52«)+e^' +5 (52'«), 

where we expressed q J+3 (c^/«) = k 2 31k and q J ' +5 (9^/c) = — (2/t<9 2 /t + K 4 )d J s K. Hence, we 
deduce 

d t d> s K = di +2 K + q J+3 (din) - 2edi +4 K - he^&^K + eq j+5 {di +1 K) , 
which gives the inductive step. □ 

Lemma 3.12. For any j G N we have 

d t f \d j s K\ 2 ds = -2 j \di +1 n\ 2 ds-Ae j \d{ +2 K\ 2 ds 

+ / <\ 2j+ \d j s K) ds + e j q 2j+6 (dl +1 K)ds. 

Proof. Using ()3.10|) . ()3.13|) and integrating by parts we deduce 
(3.14) d t j \dln\ 2 ds=2 j dlKd t d j s Kds+ j \d j s K\ 2 KE £ ds 

= 2 I & s k (dl +2 K + q j+3 (<9^)) ds 



+ e I 2&> s k ( - 2&1 +a k - hn 2 dl +2 K + q j+ \d j s +1 K)) ds 

J-y 

- / \diK\ 2 n{K-2ed 2 n-eK, 3 )ds 
= -2 [ (\dl +1 K\ 2 + q^' +4 (3>)) ds 
-As [ (\di +2 K\ 2 + q^ +e (dl +1 K)) ds. 

Proposition 3.13. For any j G N we have the e -independent estimate, fore < 1, 

// r \ 2j+3 , » \ 2j+5 

\dln\ 2 ds <C n 2 dsj +C\^jn 2 dsJ +C 

where the constant C depends only on 1/L. 



□ 
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Proof. We estimate the term f q 2j+4 (^«;) ds as in (2^1 Sect. 7]. By definition, we have 
with all the c m i less than or equal to j and 



1=1 



for every m. Hence, 



and setting 



we clearly obtain 



N„ 



|q 25 ' +4 (^)l<Eni^' 



m 1=1 



Qm = U\d:- i K\, 



1=1 



\q 2j+A (diK)\ds<J2 [ Qmds 

m T 



We now estimate any term Q m via interpolation inequalities. After collecting derivatives 
of the same order in Q m we can write 

j j 
(3.16) Q m = Y[\d l s K\ a ™ with ^« mj (z + l) = 2j+4. 

i=0 i=0 

Then 

i 



[Q m ds= /n i9>r™ ds < n c / i9^r- A ' ^ Ai = n 

T ^ ,.=0 i=0 ^ ' i=0 



S 1 1 Jj^mi 



where the values Aj are chosen as follows: Aj = if ajj = (in this case the corresponding 
term is not present in the product) and A, = - 2j (|"+i) if a mi ^ 0. Clearly, a m i\i = > 

^ > 2 and by the condition in (JUS), E^o J = E^=o £ ^f£? i = 1- 

As a m jAj > 2 these values are allowed as exponents p in inequality ()3.7|) and taking 
m — j + 1, n — i, u — k, we get 

113*11^ < ||«||^+ I ^- 7 ||«||y < C7 (\\di +1 K\\ L2 + \\K\\ L2 y- \\4\7 

with 

_ i + 1/2- l/(a mi Xi) 

J + 1 

and the constant C depends only on 1/L. 



I 1 %i 
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Multiplying together all the estimates, 
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r j 

(3.17) Q m ds< CY[(\\di +1 4 L 2 + \\4 L 2) a,, a "" ii-ii l - T 

i=0 



\ K \\L 2 



< C(\\dl +1 K\\ L2 + \\K\\ L2 ) 



^ i = Q a mi a mi II 1 1 S i — Q a rni (1 °~7ni) 



K 



L 2 



Then we compute 



^ _ ^ a mi (i + 1/2) - lj\ _ Umi{i + 1/2) - 1 

J + l J + l 



i=0 i=0 



and using again the rescaling condition in (jH.lfij) . 



2 ^ + !) 



Since 



i + 1 2 j + 4 



' 7 + 1 7 + 1 

i=0 i=0 J J 



we get 



^ ^ 2j 2 + 4j + 1 _ o 1 ^ 



8=0 



Hence, we can apply Young inequality to the product in the last term of inequality (|3.17|) . 
in order to get the exponent 2 on the first quantity, that is, 



j Q m ds < (||^ +1 k|| L 2 + ||k||l2) 2 + ||ac||^ 2 < S m J \&1 +1 k\ 2 ds + 8 m J \k\ 2 ds + ||k| 



L 2 ' 
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for arbitrarily small 5 m > and where (3 is given by 

3 1 

P = 2, a mi{^ - Vrni) —j 

i=0 1 2 



' mi 



2 n Qmi&n 

2V j a ._ 4 J+ 6 -n=o' 

Z 20+1) 
2 2(j + 1) Eio Q mt - 4j - 6 + Eio <*r 

4j+4-4j-6 + Ei= «™ 
2 (2j + 3) Ej =0 Qmi-2(2j + 3) 

2(2j + 3) . 



Therefore we conclude 



\ 2j+3 



[ Q m ds<5 m f \di +1 K\ 2 ds + S m I K 2 ds + C ( I k 2 dsj 

J 7 J 7 J 7 W 7 / 

Repeating this argument for all the Q m and choosing suitable 5 m whose sum over m is less 
than one, we conclude that there exists a constant C depending only on 1 /L and j G N 
such that 

j q 2j+4 (dlK) ds < j \d{ +l K\ 2 ds + C\^J k 2 J +C. 
Reasoning similarly for the term q 2j " +6 (£^' +1 /c), we obtain 

I q 2j+G {dln) ds< ( \dl +2 n\ 2 ds + c( ! K 2 dsY +C. 
Hence, from (j3.14j) we get 

d t j \d j s n\ 2 ds< - f \di +1 K\ 2 ds-e j \di +2 K\ 2 ds 

\2j+3 , , \2j+5 



+ C\^jK 2 dsj +Ce\^J n 2 dsj +C 

/ p \ 2j+3 /■ » \ 2j+5 

<C\^jK 2 dsj +C\^jK 2 d S J +C 



'7 / V7 

when e < 1 and the constant C depends only on 1/L. □ 
By means of Propositions 13. 101 and l3~T3*l we have then the following result. 
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Theorem 3.14. For any j G N there exists a smooth function TP ' : R — > (0, +oo) such that 

d t [ \d{K\ 2 ds <v( [ K 2 ds 



'7 Wf 

for every e < 1 and curve 7 evolving by the gradient of the functional Q £ . 

Proof. The statement clearly follows by Propositions 13.101 and 13.131 since by Lemma 13.21 
the quantity 1/L is controlled by J k 2 ds. 

The smoothness of the functions T) is obtained choosing possibly slightly larger constants 
in inequalities ()3.15|) and (|3.12|) . □ 

This proposition, like its analogue for the general case (Theorem lfi.2|) . is the key tool in 
order to get ^-independent compactness estimates. Indeed, for example, one can see that, 
by an ODE's argument, since all the flows (letting < e < 1 vary) start from a common 
initial smooth curve, fixing any j G N, there exists a common positive interval of time such 
that all the quantities ||9*k||x,2, for i G {0, . . . , j} are equibounded. This will allow us to 
get compactness and C°° convergence to the mean curvature flow as e — > 0. 

4. The General Case 

If k > [n/2] + 2 it is shown in [Hj that for every e > all the flows (p £ , associated with 
the functionals 



g' k (<p)= [ {l+e\A k \ 2 ) dfi. 

J M 



and starting from a common initial n-dimensional smooth compact immersed submanifold, 
are smooth for every positive time. 

By means of Theorem 4.5 and Theorem 5.9 in [2] and the results of ^3], the first variation 
of the functional Q e k is given by 

(k — 2)-times 

E e = — H + 2ek(-l) k - 1 (^ M oA M l ... o A A/f rY + eq 2fc - 3 (V 2fc - B B)- L 
where q 2fc_3 (B) takes values in IR n+m . 

Here we denote with A M the Laplacian of the smooth compact Riemannian manifold 
without boundary M = (M,g), where g is the metric induced on M by the immersion. 

Then we have a solution of the geometric evolution problem for any initial smooth 
immersion ipo : M — > M n+m , that is, a smooth function if 6 : M x [0, +00) — > IR n+m such 
that 

(1) the map (p e (-,t) : M — > R n+m is an immersion for every t G [0, +00); 

(2) <p E (j), 0) = (foip) for every p G M; 

(3) the following parabolic system is satisfied 

(fe — 2)— times 

dip 6 



dt 



-E £ = H + 2ek(-l) k ^A M o A m o . . . o A M h) + eq^-^V 2 * - ^) 1 . 
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5. Evolution of Geometric Quantities 

We work out some evolution equations for the geometric quantities under the flow by 
the gradient of Q\. 

In general, if a family of immersed manifolds (p(-,t) : M — > R n+m moves by d t ip = — E, 
with the field E normal, we have 

^ = 2<B,„E>, ^ = -2^(B sZ ,E)^. 

>/ 



Now for the Christoffel symbols Tfj = ( a ^ V ; J!s£.\ gis we h ave 



wtj / 9 2 e a<p\ ls / gy oe\ ls | / ay \ a?" . 

<9t \ dxidxj ' <9x/ / \ dxidxj ' cte/ / \ dxidxj ' cte; / <9t 
Then, supposing to work in normal coordinates, q x .q x . = B^- is a normal vector, hence 



dTt 3 _ I 2 E dcp\, s / dE 



Js 



01 \ OxiOxj ' 0x\ J ® y^V' q Xi I 3 

d 2 l„ d<p\ u /dE d\ \ ls (8E d\ \ , s 



Ox^Xj \ ' <9xi / \ ' OxjOxi i \ Oxj ' OxiOxi / ^ 

i / E gV \ ^ / B ; ^ \ i, 

\ ' OxiOxjOxi I \ 13 ' Ox 1 1 
= (B jh V,E) + (B a , V,E) + (E, ViBj,) - (B -, V,E> </ s . 

Remark 5.1. By this computation, since the Christoffel symbols are symmetric in the r/'- 
indices, the covariant 3-tensor (VB)- 1 is symmetric (as in the codimension one case). 

Then, we compute the evolution of B, 

dBy _0 ( gV dy? \ _ d 2 E gT^ 0^ OE 



<9t Ot \OxiOxj 13 Ox s J Ox^Xj Ot 0x s 13 0x s 

= - v « E -TC = - v i E + « VB ' VE )|l' 

where f s is the polynomial expression above in B, E and their derivatives. 
Before proceeding we need the following technical lemma. 

Lemma 5.2. If X is a vector field on M with values in M, n+m , we have 
(VtX) 1 - = + B ijg *> (X, V /0x s ) + (X, Bij) g js ^- , 

\M _ T7 vM _ R .. n js lv n, n /a„ \ _ lv R..\ n js^£_ 



(V t X) M = V,X M - Bijff" (X, Ocp/Oxs) - {X : B^) g> 



Ox, 
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More in general, 



fc-i 



(V n ... lfc X)^ =V il ... ifc X- L + X;pJ" i (V fc - i " 1 B) {WX,dtp/dx a ) 



3=0 

k-l 



Proof. We compute 

(V i X) ± = V,X X + (V.X^ - (V,X ± ) M 



= + B^ (X, - 

= V,X X + B ij9 i s (X, d<p/dx s ) + (X, Bij) g js ^ 



dx s ' 

The second formula is similar. 

The third formula follows by induction, once one works in normal coordinates where 
q x .q x . = B^, which is a normal vector. □ 

Remark 5.3. Roughly, this lemma says that the interchange of differentiation and projection 
operators introduces some extra terms in B, X and their derivatives, and the order of 
differentiation of X is lower than the initial one. 

This is useful when X is a function of B, in particular, if X is the mean curvature vector 
H we have 

fc-i 

(5.2) (V„., fc H)^ = V il ... ifc H- L + P^V^B) (V'"H, dtp/dx,) 

3=0 

+ J2(vmy-^-^B))^- 

3=0 S 

= V n ... ifc H + p^V^B) + p^V*"^)^ , 

as H is a normal vector. 

Lemma 5.4. For any s GN we have 

(5.3) |- f \V s B\ 2 dfi=-2 [ \V s+1 B\ 2 dfi+ [ q 2s+4 (V s+1 B) d/j 

&t JM JM JM 

-Ask [ |V s+fe - 1 B| 2 ^ + £ f q 2k+2s (V 2k+s - 3 B)dfi, 

JM JM 
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and 

(5.4) |- / \B\ 2s+2 d l i= [ q 2s+4 (V 2 B)dfi + e [ q 2k+2s (V 2k - 3 B) dfi . 

& J M J M JM 

Proof. Substituting E e in place of E in (|5.1|1 and expanding, after some computation using 
Lemma l5~2| we get 

(k — 2)-times 

dB ^~ 



(5.5) ^=M(-lfV 2 (A M oA M o...oA M H) + V 2 H 

+ ep^iV^B) + £ pf-i(V 2fe - 3 B)!^ + p 3 (V 2 B)|^ . 

J OXj J OXj 

Applying now formula (|5.2j) to the first term on the right hand side of (|5.5|) . we obtain 

(k — 2)-times 

dB 



^ = 2ek(-l) k V 2 A M o A A1 o ... oA M R + V 2 H 

+ ep 2k -\\7 2k - 3 B) + £ p 2k -\W 2k ~'B)^- + p 3 (V 2 B) |^ 



We now observe that for any tensor T, we have 

— VT = V— T + T * B * VE £ + T * VB * E £ 
at at 

d 

= V— T + T * p 3 (VB) + eT * p 2fc - x (V 2fc ~ 4 B) . 

L/ v 

Then, starting from equation ()5.5j) and working by induction, again in normal coordinates, 
we get 

[k — 2)— times 

d 



, V S B = 2ek(-l) k V s+2 A M o A M o . . . o A M H + V S+2 H 
at 



+ ep 2fc+s-l( V 2fc+s-3 B ) + p ^+3( VS +l B ) 

+ e pf +- 1 (V 2fc+ - 3 B)|^ + pf 3 (V S+2 B)|^ 



Hence, 



(k — 2)-times 

^|V S B| 2 =4ek(-l) k V s + 2 lJw A M oA M o ... oA'' 1IV]... ; . . . 



, 2V S+2 HV S B n iajs n lp n wz 
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Thus, we have that the time derivative of the quantity j M | V S B| 2 d\i is given by 



(k — 2)-times 



Aek(-l) k [ Vgl.fc, A M o A M o . . . o A M H V^..., B pz g** . . . g^g lp g wz 

JM 

+ 2 / VfJ^H V h „j B pz g** . . . gU-gpg™ dfi 

JM 

+ e [ q 2k+2s (V 2k+s - 3 B)dfi+ [ q 2s+4 (V s+1 B)d//, 

JM JM 

since 

—d[i = (H, E £ ) d/x, 

hence its contribution can be absorbed in the last two terms. 

We need now the following formula, which follows by direct computation, 



(5.6) ViB jt = VjBu + (B * B) 



Indeed, the 3-tensor V^B^ taking values in M n+m is not symmetric (see also Remark |5.1j) . 
Reasoning then as in Prop. 7.2-7.4], with the only care that instead of applying Propo- 
sition 2.4 in that paper, we use formula ()5.6|) . we finally obtain fl5.3j) and (|5.4|) . □ 

6. £— Independent Estimates 
For any integer s > n/2 and e > we set 

(6.1) Q S M= I (l + |V*B| 2 + |B| 2s+2 ) d/i, te[0,+oo). 

JM 

Letting e > vary, we want study the evolution of Q s e under the flows associated with the 
functionals Qt. 



(6.2) 



By the computations of the previous section we have 
dQl 



_2 I \V S+1 B\ 2 dfi - 4ek [ |V a+fc_1 B| 2 d/i 

9t J m JM 

+ / q 2s+4 (V s+1 B)rf/i + e f q 2k+2s (V 2k+s - 3 B)dfi 

JM JM 



In order to deal with the polynomial terms we state in other words Proposition 6.5 in 
(see all Section 6 in the same paper). 

Proposition 6.1. Choosing some 5 > and setting D = Vol(M) + ||H|| L n+«S(- M ), there exists 
a constant C depending only on n, m, I, z, j , p, q, r, 5 and D, such that for every manifold 
(M,g) and covariant tensor T = 7^...$,, the following inequality holds 

(6-3) ||V j T|| L p ( ^ < C \\T\\wz,v(ri\\ T \\v^) , 
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for all z G N, j G {0,...,z}, p,q,r G [1, +00) and a G \j/z,l] with the compatibility 
condition 

1 j (\ z\ I- a 
- = J - + a( + . 

p n \q n J r 

If such a condition gives a negative value for p, the inequality holds in ()6.3|) for every 
p G [1, +00) on the left hand side. 

This clearly implies, looking at the definition of the quantities Q s e , that we can alter- 
natively let the constant C in inequality ()6.3j) depend on n, m, s, I, z, j, p, q, r and 

qK2]+1 

Working now as in Section 7 of with s > n/2 fixed, we can interpolate the polynomial 
terms as follows, 

f q 2s+4 (V s+1 B) dfi < [ |V s+1 B| 2 c//i + C 1 (g[ n/2l+1 ) 

JM JM 

[ q 2k+2s (V 2k+s - 3 B)dfi<3k [ |V s+fe - 1 B| 2 rf/i + C 2 (g[ n/2l+1 ) 

JM JM 

where C^Q 1 ^ 1 ) and C 2 {Qf ,2]+1 ) are some constants depending only on n, m, fc, s and 

Q [ e /2]+1 - 

Hence, for every s > n/2, by ()6.2|) we have the estimate 

^T<- / \V S+1 B\ 2 d^i-ek [ IV^BI 2 du+C^Q^^+s^Q^ 1 ) <C{Q^ +1 )., 

& JM JM 

where C\, C 2 and C depend on e only through 

Theorem 6.2. For any integer s > n/2 there exists a smooth function Z s : R — > (0, +00) 
such that 

(6.4) d t J (1 + |V S B| 2 + |B| 2s+2 ) dfi < V (J (1 + |V [n/2]+1 B| 2 + |B| 2 ["/ 2 ]+ 4 ) d^j 

for every e G (0, 1) and any smooth evolution by the gradient of the functional Q\. 

Proof. The functions Z s can be clearly chosen to be smooth, possibly slightly enlarging the 
constants in the last inequality above. □ 

As a consequence we get the following corollary. 

Proposition 6.3. In the same hypotheses of Theorem \6.S\ there exists a continuous non- 
increasing function : (0, +00) — > (0, +oo) ; independent of e > 0, such that for every 

T eR and £ G [T, T + Q{Q [ e' 2]+1 {T))} we have Q [ e /2]+ \t) < 2Q [ e /2]+1 {T). 

Proof. The statement follows by a standard ODE argument applied to the differential 
inequality 

Q t (1 + | V [n/2]+l B |2 + | B |2[n/2]+4) ^ < z [n/2]+l H ^ + | V [»/2]+l B |2 + | B |2[n/2]+4) ^ 

which is the first case of Theorem 16.21 □ 
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In other words, this proposition says that we have an e-independent control qI™^ 2 ^ 1 < C 
in some ^-independent time interval [T, T + 0] (hence also a control the constants in 
Proposition 16. II and on the right hand side of inequalities (|fi.4j) for every s > n/2), with C 
and G depending (smoothly) only on the value of Q^^ 2 ^ +1 at the starting time T. 

7. Convergence to the Mean Curvature Flow 

In this section we prove the convergence of solutions pf : M x [0, +oo) — > R n+m to (jl.lj) 
(all starting from a common immersion <po) to the mean curvature flow ip : M X [0, T S i ng ) — > 
j^n+m De f ore its first singularity time. 

We need the following result (which can be proved as in the codimension one case as 
in [2H1 Proposition 6.3]). 

Proposition 7.1. If a manifold (M,g) satisfies Vol(M) + ||H||xn+«^ < D for some 5 > 
then for every covariant tensor S = S^...^ we have 

max |5| < C (|| VS\\ L v^) + H^Ulp^)) i/p > n, 
where the constants depend only on n, m, I, p, 5 and D. 

Let Q s (t) denote, for each nonnegative time t before the first singularity, the right hand 
side of din]) for the mean curvature flow <p at time t. 

Lemma 7.2. If the family of immersions (p £ (-,T) : M —>■ ~R n+m are bounded in the C°° 
topology, for any s G N all the quantities |V S B| are uniformly bounded by e -independent 
constants C s < +oo, in the time interval [T, T + B(sup e>0 Qe {T))], where G is the 
function in Proposition \(i., 6 A 

Proof. By the C°° boundedness of the family <//(•, T) : M — > M n+m , all the quantities 
Qe 2 ' +1 (T) are equibounded. As the function is continuous and nonincr easing, setting 
r = 6(sup e>0 Ql?^ +1 (T)) > 0, by Proposition 16 . 31 there exists a constant C > such that 
Q [ s /2]+l (t) < C for every e > and t G [T, T + r]. 

Then, again by the boundedness of the family ip e (-,T) and Theorem 16 .2\ in the same time 
interval [T, T + r] all the quantities 



/ (l + |V s B| 2 + |B| 2s+2 ) dfi 

JM 



for every s > n/2, are bounded by ^-independent constants C s < +oo. Moreover, all the 
constants in the interpolation inequalities of Propositions 16.11 and 17.11 are also bounded. 
As a first step we see that, by means of Proposition 16 .1\ we get the following estimates, 
for every p G [2, +oo) and s G N, 



|V S B| P ^ < C, 



' M 

in the same time interval [T, T+r]. Here again the constants C SiP < +oo are e-independent. 
Then, we conclude by means of Proposition 17.11 □ 
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Lemma 7.3. Assume that at time t = T the family of maps ip £ (-,T) : M — > R n+m 
converges as e — > in the C°° topology to the immersion ip? '■ M — > R n+m . Then the maps 
(f £ smoothly converge in the time interval [T,T + 6(Q [ " /2]+1 (T))) to the solution of the 
mean curvature flow starting from <pt- 

Proof. By the previous lemma, we have uniform bounds on B and its derivatives in the time 
interval \T,T + r] with r = 0(sup e>o Q^^ +1 (T)). Then, there exists C > independent 
of e such that 

dip £ {p,t) 



dt 



|E £ (p,t)| < C V(p,t) E M x [T,T + r], e > 0. 



Now we consider the metric tensors gfj (p, t) 
{v 1 } G T p M. Then we have 



dxi 



dxj 



and fix a vector V 



2|(E B ,B y )t;V| < 2|E £ | IBI^jlVl^, < C\V\^ t] 



where C does not depend on e. 

Then a simple ODE argument shows that the metrics are all equivalent; more precisely, 
there exists a positive constant C depending only on tp T such that 

(7.1) ^<^>,t)<CId, 



C 



as matrices. 

Moreover, as functions, all the g\ { = tt£- are equibounded above by a common constant. 

Hence, by Ascoli-Arzela's Theorem, up to a subsequence, the immersions tp e uniformly 
converge, as e — ► to some Lipschitz map (p : M x [T,T + r] — > ]R n+m , which clearly 
satisfies (p(p, T) = ifrip) for every p e M. 

Similarly, as the time derivative of the Christoffel symbols is given by 



(7.2) 



dt 13 



VE £ * B + E £ * VB 



(see the beginning of Section EJ) and all the metrics are equivalent, it follows that all the 
Christoffel symbols are equibounded. This means that estimating the covariant derivatives 
is equivalent to estimate the standard derivatives in coordinates, hence, we have immedi- 
ately \d s V l B\ < C s i for every s, I e N. 



Since 
we get 

and, by formula ([7.2)1 . 



0_ 

dt 



9ii 



2(E e ,B^) 







< c 
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for every s£N. 

Hence, we get |9 s ^r^-| < C s which implies, as the family of maps tp £ T is bounded in the 
C°°-topology, that |9T^| < C s . 

Since we already know that \tp £ \ are equibounded, \d(p £ \ < C and d 2 ip £ = Tdip £ + B, by the 
estimates |<9 S V'B| < C S>1 , we can conclude that the derivatives |<9 S <^ £ | are all bounded by 
e-independent constants C s , for every s G N. 

Finally, the uniform control on the time and mixed derivatives of ip £ follows using the 
evolution equation. 

Hence, the sub-convergence (p £ — > (p, as e — > 0, is in the C°° topology and (p is smooth, 
moreover, the limit metric is positive definite by (j7.1)l . 

Passing to the limit in the evolution equation dt<p> £ = E e , by the bounds on B and its 
derivatives, shows that (p : M x [T,T + r] — > IR n+m flows by mean curvature with a starting 
smooth datum tp?- Since this flow is unique, all the sequence of maps (p £ converges to (p 
which hence coincides with (p. 

Chosen now any 5 > 0, let e > be such that 

sup Q[ n/2]+1 (T) - Q^ +1 (T) < 5. 

0<e<e 

Since G is nonincreasing (see Lemma El , in the interval [T,T + Q(Q^ n W +l (T) + 5)} the 
sequence (p £ converges to the mean curvature flow (p. Letting 5 to zero, as the function 
is continuous, we get the thesis. □ 

We are now in the position to conclude the proof of Theorem 11.11 

Proof of Theorem M . 1\ Let T max be the maximal time such that ip £ converge to the solution 
of the mean curvature flow equation (p in C°°(M x [0,T max )) starting at time from the 
common immersion (po. Observe that T max is positive by Lemma \7.3\ We want to show 
that T max coincides with the first singularity time T S i ng for <p. 

Assume by contradiction that T max < T sing . Then <£>(-,£) — > <p(-,T masi ) in C°°(M) as 
t — > T max . As the function G is continuous, there exists 

lim &(Q [n/2]+1 (t)) = e(Q[™/ 2 i +1 (T max )) = e > . 

Choosing now a time T G [T max - 9/4, T max ) such that Q{Q^ +1 {T)) > 6/2, and applying 
Lemma 17.31 we see that ip £ (-,t) converges to the mean curvature flow also for t in the 
interval [T, T + 9/2}. As T + 9/2 > T max — 9/4 + 9/2 > T max , we have a contradiction. □ 
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